Anomaly in the stability limit of liquid helium 3 
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We propose that the liquid-gas spinodal line of 3 He reaches a minimum at 0.4 K. This feature is 
supported by our cavitation measurements. We also show that it is consistent with extrapolations 
of sound velocity measurements. Speedy [J. Phys. Chem. 86, 3002 (1982)] previously proposed 
this peculiar behavior for the spinodal of water and related it to a change in sign of the expansion 
coefficient a, i. e. a line of density maxima. 3 He exhibits such a line at positive pressure. We 
consider its extrapolation to negative pressure. Our discussion raises fundamental questions about 
the sign of a in a Fermi liquid along its spinodal. 

PACS numbers: 67.55.Cx, 64.60.Qb, 65.20.+W 



Because of its high purity, liquid helium is an ideal ma- 
terial in which to study homogeneous nucleation. Dur- 
ing the past years, a powerful experimental method has 
been developed to investigate the gas nucleation in the 
stretched liquid. It consists in focusing a high ampli- 
tude sound wave in a small region far from any wall, 
thus making heterogeneous nucleation unlikely. In this 
Letter, we present a detailed analysis of results recently 
obtained in 3 He M . We show that these results disagree 
to some extent with existing theoretical descriptions of 
the liquid at negative pressure. We then propose an es- 
timation of its stability limit (the spinodal line, where 
(dV/dP)T and (dV/dT)p diverge) based on sound veloc- 
ity measurements by Roach et al. || . The spinodal P s (T) 
we obtain exhibits a minimum at 0.4 K and gives a tem- 
perature dependence of the cavitation pressure consistent 
with our measurements. For water, Speedy |J previously 
proposed a spinodal with a minimum and showed that 
this change in slope of P s (T) was linked to the change in 
sign of the isobaric expansion coefficient a. This change 
in sign, which corresponds to a line of density maxima 
(LDM), also occurs in liquid 3 He at positive pressure. 
We give an extrapolation of the LDM at negative pres- 
sure and finally give theoretical arguments about the sign 
of a near the spinodal line. 



In any substance below its saturated vapor pressure, 
the liquid phase can be metastable since an energy barrier 
E\j must be overcome for liquid-gas separation to occur, 
that is for a bubble to nucleate in the liquid. The ex- 
periments in liquid helium reported in Ref. |IJ measure 
the probability of these cavitation events. For a given 
experiment performed in an experimental volume V and 
during an experimental time r, this probability is, at a 
pressure P and a temperature T: 



E(P,T) = 1 -exp 



-TqVt exp 



E h (P,T) 

knT 



(1) 



where Tq is a prefactor discussed below. The experi- 
mental measurements are reproduced by the asymmetric 
S-curve formula of Eq. ([!]) with great accuracy Q. Tq 
has the dimensions of frequency times an inverse vol- 
ume. It is natural to estimate Tq as an attempt fre- 
quency v at which the fluctuations try to overcome the 
nucleation barrier multiplied by the density of the critical 
nuclei which can be taken to be spheres of radius R c . 
Typically, R c is around 1 nm and the attempt frequency 
varies from koTjh to Eb/h; all the different estimates 
thus lie between 5 x 10 36 T and 1.5 x KpTm" 3 s" 1 K -1 . 
Pettersen et al. 0] have calculated V and t for the 
experimental method which uses an acoustic wave to 
produce a negative pressure swing in the liquid. For 
3 He and for a 1 MHz acoustic wave as in Ref. jjj , this 
gives Vt = 1.2 x 10~ 22 m 3 s ||. The theoretical esti- 
mates of the factor TqVt thus vary from 6 x 10 14 T to 
1.8 x 10 16 TK~ 1 . Although this range extends over two 
orders of magnitude, it does not significantly affect the 
value of the energy barrier: for S = 0.5, all estimates 
give E h = (34 ± 3)fc B T. 

Recently, Caupin and Balibar jjj have given experi- 
mental limits for the pressure at which the cavitation 
probability is one half (the cavitation line P ca .v(T)) in 
liquid 3 He and liquid 4 He. The upper and lower bounds 
deduced from their measurements in 3 He are shown in 
Fig. |l|(a). The actual cavitation line is located between 
these limits and parallel to them. For both isotopes, 
Caupin and Balibar made a comparison with the spin- 
odal pressure at which the nucleation barrier vanishes; 
they used the values obtained by Maris at low tempera- 
ture: P s = -9.6 bar for 4 He § and P s = -3.15 bar for 
3 He 10 1 . In fact, because of thermal or quantum 
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FIG. 1. (a) Comparison between experimental and the- 
oretical cavitation pressures. The experimental limits for 
the cavitation line are given by circles for the upper bound 
and squares for the lower one; dotted lines are guides to 
the eye. Other lines are theoretical cavitation lines calcu- 
lated with TqVt = 6 x 10 14 TK _1 using two different sources 
for the spinodal pressure (see Fig. ^): Guilleumas et al. |J 
(dashed line) and this work (solid line), (b) Temperature vari- 
ation AP = Pcav(T) — fcav(0.1 K) of the cavitation pressure. 
AP obtained from the experimental lower bound is given by 
squares; the dotted line is a guide to the eye. Other lines are 
calculated AP according to Guilleumas et al. || (dashed line) 
and this work (solid line). 

fluctuations in the liquid, the cavitation pressure is al- 
ways higher than the spinodal pressure, and the differ- 
ence can be calculated from Eq. (|]J) if the expression of 
E h (P,T) is known. 

Maris |l(| has calculated E\>{P) at low temperature 
by a density functional method; close to the spinodal, 
his results are well represented by a power law: 



fc R 



(2) 



with (3 = 47.13 Kbar~ 3 / 4 and 5 = 3/4. However, to cal- 
culate the cavitation pressure up to 0.6 K, we need to 
know the temperature dependence of Eh- The strongest 
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FIG. 2. Comparison between two theoretical estimates of 
the spinodal line: Guilleumas et al. Q (dashed line) and this 
work (full squares; the solid line is a guide to the eye). The 
spinodal found in this work shows a minimum at 0.4 K. The 
dotted line is a linear extrapolation of the LDM as measured 
by Boghosian et al. jll| between and 11 bar (see Fig. |^). 
Notice that the pressure scale is different from Fig. [lj 

source of this dependence is that the spinodal pressure 
varies with temperature; therefore we write Eb(P,T) = 
E\,{P — P S (T)) and assume that Eq. ^ remains valid at 
higher temperature with parameters (3 and 8 held con- 
stant. The temperature dependence of the cavitation 
pressure follows from Eq. (fj]): 



Pcav(T) = P S {T) 



T (T Vt 



1/5 



(3) 



Guilleumas et al. |J have calculated the spinodal pres- 
sure as a function of temperature; they find a monoton- 
ically increasing pressure as shown on Fig. |^. Guilleu- 
mas et al. have also calculated Eh(P, T) and deduced the 
cavitation line P cav (T). Their estimate of P cav (T) does 
not agree with our results. However, this estimate was 
based on a value of Vt of 2.5 x 10~ 19 m 3 s which is much 
larger than the one corresponding to our experiment. We 
have therefore used our approximation of E^ (P, T) and 
their result for P S (T) to calculate P cav (T) for the ap- 
propriate value of Vt. Fig. |l](a) shows the line obtained 
with Eq. (||) for the lowest possible value of the prefactor, 
namely T Vt = 6 x lO^TK" 1 . We first note that this 
line does not lie between the experimental limits at low 
temperature jl2|]; however, this could be due to some sys- 
tematic error in the lower bound estimate, which would 
result in shifting each pressure by the same amount (up 
to ±0.15 bar as explained in Ref. 0]). Let us focus on 
the temperature dependence of the cavitation pressure, 
which is free of this systematic error: Fig. [l](b) displays 
the quantity AP = P cav (T) - P cav (0.1K) for the three 
lower lines of Fig. 0(a). The temperature variation of AP 
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we obtain with the P S (T) curve of Guilleumas et al. is 
stronger than the experimental one. 

How can we explain this discrepancy? Of course one 
can assume that the theory fails in estimating the value of 
Tq. However, to reproduce the experimental temperature 
dependence of P cav would require TqVt to be at least 3 
orders of magnitude smaller than expected. We do not 
see any reasons to support this hypothesis. Instead, we 
think that the experimental measurements question the 
shape of the spinodal limit. 

Before proceeding further, we need to recall how 
the spinodal pressure P s can be obtained: Maris' 
method |],[l(| consists in extrapolating measurements of 
the sound velocity c at positive pressure with a law of 

1/3 

the form c = [b(P — P s )] . Maris used for c the mea- 
surements of Abraham et al. at low temperature [ jl3| . 
We used the same method with a set of data from 
Roach et al. ||: they measured the first sound velocity 
along isochores at starting pressures from 1.6 to 28.1 bar 
and as a function of temperature from 0.01 to 0.6 K. 
The spinodal line we obtained is shown in Fig. the 
spinodal pressure reaches a minimum of —2.9 bar around 
T = 0.4 K. We would like to emphasize that none of 
the previous estimates of the spinodal pressure in liquid 
3 He Hp) 14 1^] has mentioned the possible existence of 
a minimum in the spinodal line. The new shape of the 
spinodal curve we propose is sufficient to remove the dis- 
crepancy stated above: using again Eq. with the value 
ToVt = 6 x 10 14 TK _1 , we find a cavitation line which 
has a temperature dependence consistent with the exper- 
imental results (see Fig. |l|(b)). One can wonder if the use 
of Maris' method to estimate the spinodal pressure is rel- 
evant for the points at higher temperature. Thermody- 
namically, the spinodal line is the locus of points at which 
the isothermal sound velocity ct vanishes, whereas the 
measured sound velocity is the adiabatic one, cs- This 
was stated before |l6| , along with the fact that the differ- 
ence vanishes at zero temperature. We made the appro- 
priate corrections. Both spinodal lines, obtained either 
with cs or ct, show a minimum, and their difference is 
less than 30 mbar. We have actually plotted the spinodal 
obtained with the isothermal data in Fig. p], and we used 
it in the above reasoning. 

We shall now turn to the physical origin of such a min- 
imum in the spinodal. A similar behavior was first pro- 
posed by Speedy in the case of water || . This was the ba- 
sis of the stability limit conjecture introduced to explain 
anomalies of supercooled water: in his theory, the liquid- 
gas spinodal was assumed to be reentrant at tempera- 
tures below 35°C. A review of this topic also describing 
alternative theories can be found in Ref. p7|| . Following 
a thermodynamical analysis first developed in the case of 
helium Q , Speedy shows that close to the spinodal the 
sign of the isobaric thermal expansion coefficient a of the 
liquid is the same as the sign of dP s /dT. 
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FIG. 3. Temperature of density maximum (empty squares) 
and inverse of (da/dT)p around a = (full circles) 
as functions of pressure, derived from measurements by 
Boghosian et al. |ll] . The dotted line shows the extrapolation 
of the LDM used in Fig. ^ The solid line is a parabolic fit to 
l/(da/dT)p forced to vanish at the pressure of the minimum 
in the spinodal. 

Therefore, if the locus of points such that a(P, T) = 
intersects the spinodal, this results in an extremum in 
the curve P S (T). Water and 3 He have in common that 
both liquids exhibit a LDM: in some temperature range, 
they expand upon cooling. Therefore they may exhibit 
such a minimum in the spinodal. Wc have tried to adapt 
Speedy's conjecture to the case of 3 He. The measure- 
ments of Roach et al. give the expansion coefficient, but 
unfortunately they are made in a region of the phase di- 
agram where a is always negative; to obtain the LDM in 
3 He, we need to know the temperature where a vanishes. 
Therefore we used measurements by Boghosian et al. [ pd| , 
which extend to higher temperatures and agree well with 
Roach values in the region where both sets overlap; the 
result is shown in Fig. |^. A simple linear extrapolation 
of the LDM for pressures below 11 bar extends down to 
the minimum in the spinodal as shown in Fig. |^. In his 
original paper Speedy shows that the expansion co- 
efficient at the spinodal undergoes a jump from — oo to 
+00 at the temperature at which the LDM meets the 
spinodal. To find some evidence to support this predic- 
tion, we follow Speedy's analysis for water and derive 
the slope (da/dT)p around a = for each isobar in 
the measurements by Boghosian et al. This slope should 
diverge when the pressure reaches the spinodal. Fig. || 
shows that the experimental values are consistent with 
this prediction. 

We now give some theoretical arguments concerning 
the sign of a. The negative sign of a in 3 He at low tem- 
perature was first observed experimentally in 1958 by 
Lee et al. Jl9|. The same year, Brueckner and Atkins pc| ] 
pointed out how this behavior was related to the varia- 
tion of the effective mass with density. Indeed, using a 
Maxwell relation, we can write: 
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1 fdS 



a 



V \dP 



(4) 



In the Fermi liquid region, the heat capacity Cv is lin- 
ear in T and we have S = Cv = (m* /m) Cf where Cf 
is the heat capacity of the Fermi gas. Using Greywall's 
measurements of the effective mass and extrapolat- 
ing them at negative pressure as we did before |l6| ], we 
find that a given by Eq. (Q) remains negative down to 
the spinodal. Of course we should consider the correc- 
tions to the linear regime of the heat capacity and their 
evolution close to the spinodal. We see two sources of cor- 
rections. The first one is the contribution of phonons to 
the heat capacity, which varies as (T/c) 3 , where c is the 
sound velocity; this term could become important near 
the spinodal where the isothermal sound velocity van- 
ishes. However, this is relevant only for the long wave- 
length phonons: as stated by Lifshitz and Kagan [p2[ , 
the first correction to the linear dispersion gives for small 
momentum: 



uj k 2 = k 2 (c 2 + 2p\k 2 ) 



(5) 



where A is a constant. As the spinodal is approched, the 
dispersion relation thus becomes quadratic. A calcula- 
tion shows that the correction to a remains negligible at 
temperatures of interest here. We also note that, if the 
sound remains adiabatic at small k close to the spinodal, 
the use of cs, which does not vanish at P s , instead of 
ct would further reduce the phonon contribution. The 
second correction is due to the coupling of the quasipar- 
ticles to the incoherent spin fluctuations and varies as 
T 3 lnT. This effect has been studied by Greywall Q, 
who has shown that its amplitude decreases when pres- 
sure decreases; it is not clear to us if this is the case until 
the spinodal is reached, and this point requires further 
investigation. 

We will end with some remarks about 4 He. Its liquid 
phase exhibits two lines of density extrema at positive 
pressure. Clearly, it would be interesting to know how 
they extend into the metastable liquid region to deter- 
mine the shape of the spinodal line. This problem is also 
related to the behaviour of the roton minimum in the 
excitation dispersion curve and to the slope of the su- 
perfluid transition line in (T,P) coordinates at negative 
pressures. Several authors |Q] have addressed some of 
these issues, but a unified picture is still lacking. 

In this Letter, we have studied the temperature depen- 
dence of the cavitation pressure in liquid 3 He. We have 
shown that recent measurements disagree with existing 
theories. We then proposed a new picture for liquid 3 He 
at negative pressure. From the pressure and temperature 
dependence of the sound velocity in 3 He, we obtained a 
liquid-gas spinodal different from what was previously 
predicted: this new spinodal is reentrant, that is to say 
that the curve P S {T) exhibits a minimum of —2.9 bar at 
T = 0.4 K. This new feature in the phase diagram of 



liquid 3 He agrees with our measurements of the temper- 
ature dependence of the cavitation pressure. Following 
an analysis by Speedy in the case of water, we have em- 
phasized the relationship between this behavior and the 
negative expansion coefficient in 3 He. Finally we have 
given theoretical arguments to estimate this expansion 
coefficient at negative pressures. 
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